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METHOD AND APPARATUS FOR REDUCED RANK 
CHANNEL ESTIMATION IN A COMMUNICATIONS SYSTEM 



FIELD 

The present invention relates to wireless communications. More 
particularly, the present invention relates to a novel and improved method of 
reduced rank channel estimation in a communications system. 

BACKGROUND 

To improve the quality of wireless transmissions, communication 
systems often employ multiple radiating antenna elements at the transmitter to 
communicate information to a receiver. The receiver may then have one or 
more receiver antennas. Multiple antennas are desirable, as wireless 
communication systems tend to be interference-limited, and the use of multiple 
antenna elements reduces inter-symbol and co-channel interference introduced 
during modulation and transmission of radio signals, enhancing the quality of 
communications. The modeling, and thus design, of such a system, involves 
estimating several parameters of the space-time channel or link between the 
transmitter and receiver. 

The number of estimated channel parameters per transmit-receiver 
antenna pair is multiplied by the number of permutations of transmitter- 
receiver antenna pairs, creating increasingly complicated calculations and 
decreasing estimation quality. Therefore, it is desirable to have methods of 
channel estimation that use a reduced set of parameters. Similarly, there is a 
need for an improved method of channel estimation for radio communications 
systems having multiple transmitter antennas. 
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SUMMARY 

The presently disclosed embodiments are directed to a novel and 
improved method and apparatus for estimating channel parameters in a 
5 communication link in a wireless communication system having multiple 
transmitter antennas using a reduced rank estimation method. Each path from 
a transmitter antenna to the receiver constitutes a channel within the link. The 
number of channels, therefore, increases with the numbers of transmitter 
antennas and receiver antennas. The method exploits redundant and /or a priori 

10 knowledge within a system to simplify the channel model used as a basis for 
the estimation calculations and to improve the estimation quality. In one 
embodiment, a covariance matrix is calculated and analyzed to determine if the 
number of channel parameters may be reduced for channel estimation. If not, 
all parameters are estimated, otherwise a reduced rank channel model is used 

15 for the calculation of channel parameter estimates. 

In one aspect, a method for modeling a link in a wireless communication 
system, the system having a transmitter having N antennas and a receiver 
having M antennas, each path from one of the N transmitter antennas to the M 
receiver antennas comprising a channel, includes determining a matrix 

20 describing parametric relations of the link; ranking the matrix; determining if 
the rank is less than NxM; if the rank is less then NxM performing an extraction 
of a subspace of the matrix; deriving channel impulse responses for each 
channel based on the extracted subspace of the matrix; and demodulating a 
received signal using the channel impulse responses. The matrix may be a 

25 covariance matrix describing the link, wherein the covariance matrix represents 
a plurality of impulse responses between the transmitter and the receiver. 
Alternatively, the matrix may be a sample matrix describing the link. 

Further, determining the matrix may include estimating a plurality of 
parameters describing at least one channel. The parameters may include a 

30 distance between transmitter antennas. In one embodiment, the parameters 
include a transmittal angle with respect to a configuration of the transmitter 
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antennas. In an alternate embodiment, the matrix describes parametric 
relations of the link in the frequency domain. 

Further, ranking the matrix may include determining an eigenvalue for 
the matrix. In one embodiment, if the rank is equal to (NxM) a set of correlated 
5 impulse responses is applied for demodulating. In one aspect, a wireless 
apparatus is operative to model a link in a wireless communication system by 
determining a matrix describing parametric relations of the link; ranking the 
matrix; determining if the rank is less than NxM; if the rank is less then NxM 
performing an extraction of a subspace of the matrix; deriving channel impulse 

10 responses for each channel based on the extracted subspace of the matrix; and 
demodulating a received signal using the channel impulse responses. 

In another embodiment, a wireless communication apparatus includes a 
correlator operative to estimate a covariance matrix representing a link with a 
transmitter based on signals received from the transmitter; a rank analysis unit 

15 coupled to the correlator and operative to estimate a rank of the covariance 
matrix; and a channel estimation unit coupled to the rank analysis unit and 
operative to generate a reduced rank channel estimate. The covariance matrix 
may represent a plurality of impulse responses between the apparatus and the 
transmitter. In one embodiment, the rank analysis unit is operative to 

20 determine an eigenvalue corresponding to the covariance matrix and is 
operative to compare the estimated rank of the covariance matrix to a 
predetermined full value. 

In still another embodiment, a method for estimating a link in a wireless 
communication system includes estimating a covariance matrix for the link; 

25 determining if the rank of the covariance matrix is reducible; reducing the rank 
of the covariance matrix; and estimating a set of impulse responses for the link 
using the reduced rank covariance matrix. Additionally, the method may 
include determining a correlation of the channel; ranking the covariance matrix; 
and performing an extraction of a reduced rank matrix out of the covariance 

30 matrix. 
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In one embodiment, a wireless communication apparatus is operative 
within a wireless communication system having a transmitter having N 
antennas and a receiver having M antennas, each path from one of the N 
transmitter antennas to the M receiver antennas comprising a channel. The 
5 apparatus includes a first set of computer readable instructions operative to 
determine a co variance matrix describing the link; a second set of computer 
readable instructions operative to rank the covariance matrix; a third set of 
computer readable instructions operative to determine if the rank is less than 
NxM; a fourth set of computer readable instructions operative to perform an 

10 extraction of a reduced rank matrix out of the covariance matrix if the rank is 
less then NxM; a fifth set of computer readable instructions operative to derive 
channel impulse responses for each channel based on the reduced rank 
covariance matrix; a sixth set of computer readable instructions operative to 
demodulate a received signal using the channel impulse responses. The 

15 apparatus may further include an equalizer operative in response to the sixth 
set of computer readable instructions, wherein a configuration of the equalizer 
is determined by the rank of the covariance matrix. In one embodiment, the 
apparatus includes a seventh set of computer readable instructions operative to 
derive a correlated channel impulse response. 

20 In still another aspect, a wireless communication apparatus includes a 

channel estimation means operative to estimate a covariance matrix 
representing a link with a transmitter based on signals received from the 
transmitter; a rank analysis unit coupled to the correlator and operative to 
estimate the rank of the covariance matrix; and a channel estimation means 

25 coupled to the rank analysis unit and operative to generate a reduced rank 
channel estimate. 

Further in another aspect, a wireless communication apparatus includes 
a correlator operative to estimate a covariance matrix representing a link with a 
transmitter based on signals received from the transmitter; a rank analysis unit 
30 coupled to the correlator and operative to estimate the rank of the covariance 
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matrix; and a channel estimation means coupled to the rank analysis unit and 
operative to generate a reduced rank channel estimate. 

In yet another aspect, a method for estimating a link in a wireless 
communication system includes estimating a covariance matrix for the link; 
5 determining if the rank of the covariance matrix is reducible; reducing the rank 
of the covariance matrix; and estimating a set of impulse responses for the link 
using the reduced rank covariance matrix. The method may further include 
determining a correlation of the channel; ranking the covariance matrix; and 
performing an extraction of a reduced rank matrix out of the covariance matrix. 

10 In another embodiment, a wireless apparatus include channel estimation 

means operative to determine significant delays and determine a set of 
estimates of full dimension channel parameters associated with the significant 
delays, wherein each one of the set of estimates corresponds to an instance in 
time; eigenvalue computation means operative to determine eigenvalues of the 

15 set of estimates of the full dimension channel parameters and find any 
dominant eigenvalues; and channel estimation means operative to determine a 
set of reduced rank channel parameter estimates in response to the dominant 
eigenvalues. Further, the apparatus may include eigenvector computation 
means operative to determine at least one eigenvector associated with one of 

20 the dominant eigenvalues of the set of estimates; wherein the channel 
estimation means uses the at least one eigenvector to project the set of estimates 
of the full dimension channel parameters onto the subspace spanned by the at 
least one eigenvector. 

25 BRIEF DESCRIPTION OF THE DRAWINGS 

The features, objects, and advantages of the present invention will 
become more apparent from the detailed description set forth below when 
taken in conjunction with the drawings in which like reference characters 
30 identify correspondingly throughout and wherein: 
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FIG. 1 illustrates configurations of wireless communication 
systems including multiple transmitter antennas; 

FIG. 2 illustrates a model of a wireless communication system 
according to one embodiment; 
5 FIG. 3 illustrates a model of a channel between transmitter and 

receiver in a wireless communication system; 

FIG. 4 illustrates the physical layout of antennas in a transmitter of 
a wireless communication system; 

FIG. 5 illustrates a flow diagram of a method of reduced rank 
10 channel estimation for a wireless communication system according to one 
embodiment; 

FIG. 6 illustrates a plot of the estimation gain of one embodiment; 
FIG. 7 illustrates a system configuration according to one 
embodiment. 

15 and 

FIG. 8 illustrates an exemplary embodiment of a wireless communication 
system. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

20 

Multiple radiating antennas may be used to improve transmission 
quality in a wireless communications system. In the design of third generation 
mobile radio systems, for example, various transmitter antenna diversity 
techniques are presented. Multiple transmitter antennas may be used to 

25 communicate information to a receiver using a single or multiple receiver 
antenna (s). Multiple antenna systems offer an improvement in quality. 
However, the improvement is dependent on the accuracy of the channel model 
used in the receiver to demodulate the transmitted information. Modeling of 
the transmission channel uses parameter estimates and determines an effective 

30 channel impulse response for the channel. When multiple antennas are used, 
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the modeling involves estimates of each transmission channel for all 
transmitter-receiver antenna pairs. 

The transmission channel from transmitter to receiver is a space-time 
channel described generally by at least one impulse response. Often there is 
5 little change in the channel parameters from one channel to another, such as 
where the channel impulse responses differ only in phase. In such a case, it 
may not be necessary to derive estimates of impulse responses independently 
for each channel, but rather some information may be reused. When channels 
are correlated, a reduced rank representation of the channels may be used. 

10 Reduced rank refers to the reduced number of completely uncorrelated 
channels used to describe the link between transmitter and receiver. One way 
to observe this reduced rank is the rank reduction of the channel covariance 
matrix used to describe the mutual statistical dependencies of the different 
channel impulse responses. Note that the reduced rank can also be realized by 

15 other parameter measures. For example, in one embodiment a sample matrix is 
formed of columns comprising samples of channel impulse response estimates 
over time, wherein the reduced row rank of such a sample matrix is applied as 
described herein. A reduction in rank may result in a less complex filter or 
demodulator, i.e., reduces the number of filters and/or filter elements and/or 

20 demodulation units used in the receiver. Furthermore, the reduction of the 
number of estimated parameters used to characterize the channel leads to 
improved accuracy of the channel model. 

FIG. 1 illustrates configurations for wireless communication systems 
having multiple transmitter Tx antennas. Two paths are illustrated: a first 

25 multiple input, multiple output (MIMO) and a second path multiple input 
single output (MISO). The MISO configuration places multiple Tx antennas in 
communication with a single Rx antenna. The MIMO configuration extends 
this to multiple Rx antennas. The channel model for one of the systems of FIG. 
1 is illustrated in FIG. 2, in accordance with one embodiment, specifically, for a 

30 wireless system employing coherent demodulation, having a link between a 
transmitter and a receiver, wherein at least the transmitter employs multiple 
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antennas. The wireless communication system 10 includes a transmitter 12 and 
receiver 16 that communicate via an air interface. A channel model 14 
represents the channels for antenna pairs between transmitter 12 and receiver 
16. Channel model 14 considers the channels within a link, such as the MISO 
5 link of FIG. 1. 

Continuing with FIG. 2, let N Tx be the number of antennas used at the 
transmitter 12 and A^ the number at the receiver 16, respectively. In general, 
for each significant propagation delay between transmitter and receiver, 
( N Tx • ) transmission channels exist for the pair, wherein for a significant 

10 propagation delay the received signals resemble the known transmitted signals 
with high certainty. In other words, define N E as the number of significant 

propagation delays, also referred to as echoes. The ( N Tx ■ • N E ) channel 

impulse response samples are then estimated to perform coherent 
demodulation. When the channels are uncorrected, the ( A^ - A^ ■ N E ) channel 

15 impulse response samples are modeled as completely uncorr elated random 
processes and the estimates of these channel impulse response samples may be 
derived independently without loss of demodulation performance. However, if 
the ( N Tx - • N E ) channel impulse response samples are not uncorrected 

random processes the ( N Tx • • N E ) channel impulse response samples may 

20 be modeled as a linear combination of a smaller number N Ch of channel impulse 

response samples, wherein N Ch < (N Tx • N E ) . Such cases include, but are 

not limited to, minimal angular spread at the transmitter and/or receiver in the 
effective channels due to propagation conditions. If N Ch is known, or 

estimated, and the linear transformation of the N Ch channel impulse response 

25 samples are resolved into the ( N Tx • ♦ N E ) channel impulse response 

samples, then modeling may be accomplished with the N Ch channel impulse 

response sample estimations. This reduces the number of parameters to be 
estimated while increasing estimation quality, yielding an increase in the 
demodulation performance. Even if the exact representation of the linear 
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transformation of the N Ch channel impulse response samples into the 
corresponding ( N Tx • • N E ) channel impulse response samples is not known, 
the modeling may still be accomplished with N Ch channel impulse response 

sample estimations if the subspace spanned by vectors of this linear 
5 transformation is known or can be estimated. 

This principle is referred to as "Reduced Rank Channel Estimation." The 
transformation of the N Ch uncorrected channel impulse responses into the 

( N Tx • ■ N E ) correlated channel impulse responses can depend on factors 

including, but not limited to, the antenna configuration, antenna patterns, 
10 polarization characteristics, propagation conditions and more. In some cases 
the transformation might be known a priori, in other cases it can be derived or 
estimated, for example by angle-of-arrival estimation. The subspace spanned 
by the linear transformation of the N Ch channel impulse response samples into 

the corresponding ( N Tx - - N E ) channel impulse response samples can be 

15 determined by estimating the rank and the eigenvectors of the ( N Tx • - N E )- 

dimensional covariance matrix of the ( N Tx • • N E ) channel impulse response 

samples. This subspace can also be determined by using a singular value 
decomposition of a matrix holding columns with all ( N Tx • N E ) channel 

impulse response sample estimates for different points in time. Note that if the 
20 channel impulse response samples are corrupted by a known correlation noise, 
and if the noise correlation can be estimated, the ( N Tx • • N E ) channel 

impulse response samples may be filtered by a noise de-correlation filter first. 

In one embodiment, the rank reducing transformation is known a priori 
or is estimated. In other words, the mapping of the N Ch channels onto the 

25 {N^'N^) channels is ascertainable. The reduced rank channel is then 

estimated using the ascertained transformation. When desired, an equivalent 
full dimensional channel model may then be derived from the reduced rank 
estimate by transforming the reduced rank estimate back to the larger 
dimension. 
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In an alternate embodiment, the rank reducing transformation is not 
directly known, but the subspace spanned by the transformation may be 
extracted from the dominant eigenvectors of the channel covariance matrix. 
Note that the subspace may be referred to as the signal subspace or the channel 
5 subspace. The process involves first estimating a channel covariance matrix 
and finding the dominant eigenvalues. By determining the associated 
eigenvectors which span the channel subspace, the process projects the 
conventional channel estimate into the channel subspace, yielding a reduced 
rank channel model with reduced estimation errors. If desired, the reduced 
10 rank model may be transformed back into an equivalent full dimension channel 
model. 

FIG. 3 illustrates a model 18 of a MIMO channel for continuous time 
having a linear MIMO filter 20 with N Tx inputs and outputs. The linear 

MIMO filter 20 is defined by the x matrix H (t) comprising of linear 
15 functions h^(t) 9 i = l...N Tx , j = L.-AT^ . Generally, * ~ l.-.N Tx , j = L.JV^ are 

unknown linear functions. The linear MIMO filter 20 represents the (N Tx • N^) 
radio channels through which the N Tx transmit signals pass to the receiver 
antennas. These radio channels are characterized by their channel impulse 
responses hy(i) 9 i = 1.. .N Tx , j = 1..^^. The input signal to the model, Jc(0, is a 

20 ( N Tx x 1 ) column vector representing the N Tx band-limited transmit signals, and 
the output signal from the model, $(t) , is a ( xl )column vector , sampled at 
t = 7\ 27\.., as illustrated by switch T, where the bandwidth of the transmitted 
signals is less or equal to l/T . The received signals contain additive 
perturbation signals represented by the x 1 column vector %(t) , introduced 

25 due to noise or co-channel interference. The additive perturbation signals are 
added at summation nodes 22. The relation between the input signals Jc(/) , the 

channels H(t) , the perturbation %(t) and the output signals ^(t) is given by 

fto = tf T (f)*fto + ftr), (1) 

30 where * denotes the convolution. 
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FIG. 4 illustrates the physical configuration of antennas at the transmitter 
of an exemplary embodiment modeled as in FIG. 2. A reduced rank method is 
applied to estimate the link represented by channel model 14, having a 
transmitter 12 configured with the four (4) antennas, each spaced at a distance 
5 "d." The specifics of the configuration and model are discussed hereinbelow. 
Note that the estimation procedure is performed at the receiver 16. A reference 
direction is given by the horizontal line. Angles of transmission are measured 
with respect to this reference. The angle "a" corresponds to an angle of a 
propagation path with respect to the reference within a 2-D plane as illustrated. 

10 A range of angles with respect to the reference is also illustrated. The following 
method is used at the receiver 12 in system 10 to estimate the link. 

FIG. 5 illustrates a flow diagram of an exemplary method of channel 
estimation used to process signals in a receiver unit in accordance with one 
embodiment. Process flow begins by searching for significant propagation 

15 delays in the channel, i.e. searching for significant echoes at step 40. In one 
embodiment the process involves a sliding correlation of the received signals 
with known transmitted signals or known components of the transmitted 
signals. Correlation refers to the degree with which the received signals are 
related to the known transmitted signals, wherein a perfect correlation proves a 

20 relationship between the signals with high confidence. For time-shifted signals, 
wherein sliding delays are used to shift the received signals in time, a resultant 
sliding correlation provides the degree of certainty with which the time-shifted 
signals resemble the known transmitted signals. Thus in the wireless system 
context, sliding correlation relates to the synchronization of known signals 

25 transmitted by the Tx antennas with time-shifted versions of the received 
signals. The exemplary embodiment of reduced rank channel estimation uses 
sliding correlation of the received signals with known transmitted signals to 
estimate the number N E and the values r l9 r 29 K ,t Ne of significant propagation 

delays, i.e. delays for which the received signals shifted back by these delays in 
30 time resemble the known transmitted signals with high certainty. The 
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procedure of sliding correlation in order to find significant propagation delays 
is also known as "searching" in CDMA systems. 

The method then estimates parameters for multiple observable channels 
between the N Tx transmitter antennas and the receiver antennas at step 42. 

5 The channels are radio network connection pairs coupling at least a portion of 
the N Tx transmitter antennas to at least a portion of the receiver antennas. 
In the exemplary embodiment, there is a connection between each transmitter 
12 antenna and each receiver 16 antenna, resulting in ( N Tx • ) channels. The 

parameters describing the multiple channels are those characteristics that 
10 impact the impulse responses of the channels. Assuming that N E significant 
propagation delays (echoes) exist between transmitter and receiver, 
( N Tx • • N E ) complex samples of the ( N Tx ■ ) channel impulse responses 

could be used as a set of parameters describing the multiple channels. This set 

of parameters is denoted by a ((N Tx • • N E )xl) vector termed £ herein. The 

15 relation between %(t) , )t , ^),and $(t) is developed hereinbelow. 

With rj,r 2 ,K,r^ being the significant propagation delays between 
transmitter and receiver, the model described by (1) can be expressed as 

ft/) = f;# T (r < )-£(r-r,) + ?(0. (2) 

e=l 

This can be transformed into 

20 

fto = [(i ( ^>®^ (i^-^^a-r^^K+fto, (3) 

where ® denotes the Kronecker tensor product, 1 (Nrx) is a (N^xN^) identity 
matrix, and the vector K is obtained from the matrix H(t) such that K 
25 K total (r e ) = vect{H(r e )}, e = IK N E , (5) 

hold. The ((N Tx • N^xl) vector % totaJ (T e )is consisting of the elements of the 
matrix H(t) sampled at r e with all columns of H(r e ) stacked on top of each 
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other in the vector h total (r e ) , which is denoted by the operator vect{H{r e )} in (5), 
Ko,a,(?e) is given by 



It 

It 



M 

(«■. ) *W (T e ) K Wt (T e )]' 



(6) 



Since the output signals are sampled at a sampling rate of 1/T , vectors 

containing the discrete time samples can represent segments of a finite duration 
of the continuous time signals. For the sake of simplicity, the received signals 
$(t) are described herein by a discrete time representation over a finite 

duration of time t = 0,T,K , (Af r -l)r, where N T is the number of samples 

taken over time. Therefore, the following abbreviations are used. Each discrete 
time transmitted signal at antenna n delayed by r is given by a vector 

$ (n \r) = [x n (0-T)x n (T-r)K x n «N T -l)T-T)] T , n = lKN Tx . (7) 

Wherein the matrix describing all of the discrete time transmitted signals 
delayed by r is given as 

S(r) = [P 0) (r) F (2) (t) K P^'V)] . (8) 

The matrix A describes all discrete time transmit signals having significant 
delays. 

A = [(l (Wb) ® ^r,))^^ ®S(r 2 ))K ® S(r^))] (9) 
A vector describing each discrete time perturbation signal at antenna n is given 



as 



#»> = [ Zn (0)z n (T)K z n ((N T -l)T))\ 



n = lKN Rx . (10) 



and the vector of all the discrete time perturbation signals is given as 

K = [£ (1)T £ (2)T k £<"«* )T ] T . 



(ii) 



The vector of the discrete time received signal at antenna n is given as 

¥ W = [y„(0) y n (T)K y n «N T -1)T)] T , n = lK N*. (12) 

25 and the vector of all discrete time received signals is given as 
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Y = [^>Tp( 2 )T K ^)T]T (13) 

Using the above abbreviations, the discrete time output signals of the MIMO 
channel model 18 illustrated in FIG. 3 over a period of time from 
t = 0, T, K , (N T - l)T may be reduced to the simple model 

5 P=A-£+£. (14) 

The second step in the flow diagram in FIG. 5 at step 42, is to repeatedly process 
estimates for a set of parameters characterizing the multiple channels between 
transmitter and receiver. For the above-described mathematical representation 
10 of the channel model, this may be equivalent to processing estimates 

n = lKN h of the vector It in (14) for N h different points in time. A 

conventional method uses the correlation of the received signals, shifted back in 
time by certain delays, with known transmitted signals, such as pilot signals 
specific to the transmitter antennas, or predetermined training sequences. As 
15 the significant propagation delays t 19 t 2 ,K ,t Ne are already determined in step 

40, the exemplary embodiment of reduced rank channel estimation uses the 
correlation of known transmitted signals with versions of the received signals, 
shifted back in time by r,,r 2 ,K 9 t Ne , to generate a channel model, such as 

channel model 14 of FIG. 2, characterized by the vector ft . If the noise vector K 
20 represents spatial and temporal white perturbation, wherein the noise 

covariance matrix is given as R n = ($? H ^ = a 2 -\ {Nrx Nt) , and if the matrix A 

comprises of the a priori known signals, such as pilot symbols of a CDMA 
system, channel estimates obtained by correlation can be described by 

li = A H • Y . (15) 
25 If the noise vector K does not represent spatial and temporal white 

perturbation, the channel estimates obtained by correlation can be described by 

f=A H /?; 1 .P=A H /?; , A.f+A H /?;^. (16) 

Note that R n might be known a priori or could be estimated from the received 
signals. The channel estimate of (16) contains a perturbation vector A H R~ l ii 
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with a covariance matrix of R p = A u R~ l A. This covariance matrix is not 
diagonal in general, i.e., the components of the perturbation vector contained in 

jf are correlated in general. If R p is known or can be estimated, the 

components of the perturbation vector contained in fi could be de-correlated 
5 by transforming fi with R~ ul . This will be assumed in what follows, wherein 



%=r: u2 a h r: 1 ? 



(17) 



shall hold. 

As illustrated in FIG. 5, a covariance matrix of the channel parameters is 
estimated at step 44. Covariance measures the variance of one random variable 
10 with respect to another. In this case, the covariance matrix describes the 
variance of the various channel parameters with respect to each other. 
According to the above-described mathematical representation of the channel 

model, step 44 corresponds to processing an estimate R h of the channel 
covariance matrix R h = (fi fi H ^ . Such an estimate may be given as 

15 R^-L^fa-fa" . (18) 

If the MIMO channel has a reduced rank wherein N Ch < (N Tx • -N E ), i.e. the 
(N Tx • * N E ) MIMO channel impulse response samples can be described as a 
linear combination of N Ch uncorrected channel impulse response samples. The 

channel vector )rt can be modeled as a linear transformation of a channel vector 
20 £ of reduced dimension, wherein 

%=B-£, (19) 
and wherein B is a {{N Tx - N ^ • N E )x N Ch ) matrix describing the linear 

transformation. As given hereinabove, the vector £ is a ( N Ch xl) vector with 
uncorrelated components, i.e. R g =(F F H ) * s a diagonal ( N Ch x>N Ch ) matrix. In 

25 this case, the channel covariance matrix is given as 

R h =B R g B H . (20) 
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As a consequence, the rank of the channel covariance matrix R h is equal to N Ch . 
Given (20), and assuming that correlation according to (17) is used to derive the 

channel impulse response estimates t, , the covariance matrix of fi is given as 

R h =(%'% H ) = R l p nH B-R g -B H R p ,2 +V N * (21) 



Due to the reduced rank N Ch of R h , the eigenvalue decomposition 

R l p /2H BR g -B H R p 12 =R 1 p /2H R h R 1 p ' 2 =E'A-E H , (22) 

10 yields only N ch non-zero eigenvalues, where A is a diagonal matrix containing 
the eigenvalues and £ is a square matrix containing the eigenvectors of 
R l ' m B- R e • B H R l J 2 . With (21) and (22) the covariance matrix estimate R h may 
be expressed by 

R. = E • (A + l iN *' n **' Ne) ) • E H , (23) 



15 i.e., R h shares the eigenvectors with R p B • R g • B R p . Since A is a diagonal 
matrix with only Af c/I non-zero elements, (iV rjc • -N E )- N Ch eigenvalues of 

7? ft are constant, and A^ c/I eigenvalues of R h are larger than the former ones. 

These larger eigenvalues are termed dominant eigenvalues in the sequel. With 
a diagonal matrix A c , containing all dominant eigenvalues of the estimated 

20 channel covariance matrix, the matrix E c , containing the corresponding 

eigenvectors, and with the matrix E N , containing the remaining eigenvectors, 

(23) becomes 

R h =E c -A c -E } *+E N E%. (24) 
Therefore, the matrix E c contains the eigenvectors spanning the channel or 
25 signal subspace. 

The estimated covariance matrix R h is then ranked at step 46, meaning 
that the number of dominant eigenvalues is estimated. The rank is compared to 
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a maximum value "MAX" at step 48. MAX is equal to the total number of 

estimated channel parameters in the vector fi . In other words, MAX is equal to 
(N Tx • • N E ) . As many of the mechanisms impacting correlation, such as 

the directionality of the propagation paths, do not change quickly over time, the 
5 correlation characteristics may be estimated by averaging over rather long time 
intervals in comparison to the inverse fading rate of the channel(s). 

The rank of the covariance matrix determines whether the 
(N Tx ■ • N E ) channel parameters describing the (N Tx - ) existing 

transmission channels can be modeled as a linear combination of a smaller 
10 number N ch of equivalent uncorrected channel parameters If a reduced rank 

is available, the channel subspace E c of the estimated covariance matrix R h is 
derived at step 52. Note that instead of using the estimated covariance matrix 
R h , the rank of R h and the channel subspace E c can also be derived from the 
matrix of channel parameter estimates 

15 X h =p l) £ (2) K (25) 

by using singular value decomposition. 

With the channel subspace E c , reduced dimension channel parameter 

vectors are estimated at step 54, according to 



20 f (n) = E" • £ (n) , (26) 

effectively projecting the originally estimated channel parameters into the 
channel subspace. This projection into the channel subspace reduces the 
estimation error. If a reduced complexity demodulator is used in the receiver, 
which uses the reduced rank channel, i.e., takes only a reduced number of 
25 channel parameters into account for demodulation, the estimates of (26) may be 
directly used in the demodulator for coherent demodulation. In other words, 
processing would flow directly from step 54 to step 58, or at a minimum step 58 
would used the reduced rank estimates. 
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If a conventional receiver, designed for the full rank channel model, is to 
be used, the estimates § {n) may be transformed back into the full dimensional 
space at step 56, according to 

kl=R-;' 2H Ec-% (n) , (27) 
5 wherein the factor R~ U2H is used to make the estimate unbiased. Note that the 
estimate of the channel subspace E c may be updated continuously by using a 

sliding time window for the estimates R h or X h , respectively. This eliminates 

the delay of waiting for a new complete sample set, by using a portion of the 
previous sample set with incrementally time-shifted new values. 
10 If rank reduction is not possible, processing continues to use the full 

rank of the system to model the channel at step 50. In this case the method 
estimates the (N Tx • • N E ) channel parameters independently from each 

other. Once the system is modeled, signal demodulation continues at step 58. 

The MISO path illustrated in FIG. 1 is provided as an exemplary 
15 embodiment. As illustrated, the transmitter, Tx, has four (4) radiating antennas 
(N Tx =4) and the receiver, Rx, has one (1) antenna (N^ =1). For simplicity, 

several assumptions allow a straightforward analysis demonstrating the 
applicability of the exemplary embodiment to modeling a system as illustrated 
in FIG. 1. First, the example assumes that each Tx antenna transmits a pilot 

20 signal specific to that antenna, wherein the antenna-specific pilot signal is time- 
aligned and orthogonal to the pilot signals of the other Tx antennas. 

Second, assume the channels are frequency non-selective fading 
channels, each made up of a large number, P, of radio network paths. The 
paths each have approximately a same rim length and a same attenuation. The 

25 second assumption ensures that the relative propagation delay is smaller than 
the inverse of the transmission bandwidth. The propagation delay of two radio 
paths is typically due to differences in run length. 

Third, the channel model is restricted to 2-D propagation, i.e. all effective 
radio paths are located in a 2-D plane. See FIG. 4. Additionally, the geometry 

30 of the effective radio paths at the transmitter is assumed to be time-invariant, 
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wherein each path departure angle, measured with respect to a reference 
direction of Tx, are concentrated around an average angle, a . The radio path 
angles are Gaussian distributed having mean a and a standard deviation a . 
For one simulation, a is selected randomly between -60 and +60 degrees. The 
5 standard deviation a is assumed to be square root of two degrees. Fourth, the 
arrival paths at Rx are assumed uniformly distributed between 0 and 360 
degrees to consider local scattering. Fifth, no line of sight exists. 

Sixth, assume a specific phase and Doppler shift for each path. The path- 
specific phase is selected randomly according to a uniform distribution between 

10 0 and In . Additionally, the path-specific phase is adjusted for each Tx antenna 
according to the geometrical antenna configuration, i.e., the antenna location 
with respect to a reference point. For phase adjustment, assume object 
scattering is considered in the far field. The channel-specific Doppler shifts are 
generated according to a uniform distribution of the angles of arrival paths at 

15 Rx, a carrier frequency and a predetermined Rx speed. In the exemplary 
embodiment, the carrier frequency is assumed to be 1.8 GHz and the receiver 
speed equal to 60 km/h yielding a maximum Doppler shift of 100 Hz. In the 
exemplary embodiment, each Tx antenna covers a 120-degree sector, with the 
antenna patterns all oriented towards a = 0 . 

20 Given the exemplary system as detailed, application of the process of 

FIG. 5 provides a channel model having a time variance according to the classic 
Doppler spectrum. It is possible to consider an antenna-specific radiation 
pattern. With this channel model, the channel impulse responses for the 
channels seen through the different transmitter antennas can be generated 

25 using the same set of radio paths, thus, introducing realistic correlation in the 
fading of the different channels. 

On the receiving side of the air-interface, at the single antenna of Rx, the 
method derives an impulse response estimate for each of the four transmission 
channels, i.e., the four radio network connections between Tx antennas and the 

30 Rx antenna. The estimate is based on the a-priori knowledge of spreading codes 
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used to generate the antenna-specific pilot signals associated with each Tx 
antenna. 

Referring again to FIG. 4, in the geographical configuration of antennas 
at Tx., the antennas are positioned in a line having constant spacing d between 
neighboring antennas, wherein d = X , i.e., antennas are spaced one wavelength 
apart. Note that Rx has a single omnidirectional antenna. A total number of 
effective radio paths is considered with P = 50 . Channel specific variables, a p , 

f p , and O p , represent, respectively, the angle measured from the reference line, 

the Doppler shift and the phase. The equation describing the channel impulse 
response for a channel between Tx antenna n and the Rx antenna is given as 



f 



J * 



v 



V 



JJ 



(28) 



wherein g n (a) is the antenna-specific complex azimuth radiation pattern of each 



Tx antenna. 



If a channel has no angular spread, and all path-specific angles a p are 
15 equal to a , the channel impulse response for each Tx antenna is given as 



f (w — -)d ^ 1 p i / \\ 
h n {t) = g n (a) exp j ± — ^sin(a) •-r=S ex PU*( (I> /> +2*tf p t)yS{T-T 0 ) 

A ) ^p p=l 

in-i)d , 0| 



V 
r 



(29) 



J 



where h(t) is the equivalent channel impulse response for a equivalent isotropic 

20 Tx antenna at the reference point. 

In this case the channel impulse responses for the different Tx antennas 
only differ by a complex factor, i.e., the channels are completely correlated. The 
steering vector is then defined as 
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8{a) = 



g,(a)-exp 



( . -3d . ,_ X Y| 



J- 



v 



2A 



sin 



[a) 



J 



g 2 (tf)exp 



v 



J 



V 



j—un(a) 
d . 

— rsinfa 
g 4 (a)-exp| >~sin(Zr) 



g 3 (a)exp 



V 



(30) 



and the channel impulse response vector as 



f(r)« 



^(0 

mo. 



(31) 



The four (4) channel impulse responses seen from the Tx antennas are 
then copies of the channel impulse response hit) , weighted by four different 

complex factors, which means, the vector ti(t) is a linear transformation of the 
10 scalar h(t) given by 

%(t) = 8{a)-h(t), (32) 
i.e., the vector £ in the linear transformation of (19) is in this example equal to 

the scalar h(t) and the matrix B is equal to the vector £(a). This means the 

15 channel covariance matrix R h =(fi-ti H } is equal to R h = S(a)S(a) R (^h(t)\ 2 ^ in 

this example. If the steering vector Kific) is known, such as a-priori knowledge 
of the antenna configuration and the radio path direction a , it is sufficient to 
estimate the scalar h(t) and either calculate an estimate for using the linear 
transformation with or use the estimate of hit) and &(a) directly for 

20 demodulation. 

Note that for the case when£(or) is known it may be sufficient to 

estimate h{t) and then compute an estimate of n(t) from the scalar estimate of 
h(t) . If the demodulator is designed such that the channel consists of a single 
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scalar, i.e., the demodulation considers &(a), then it is possible to demodulate 
using K(a) and the scalar channel. 

The antenna-specific pilot signals at the transmitter are termed x n (t), 
n = IK N Tx , and the relationship is defined by 



\ Xn (tf=l V ne{lKN Tx } 



(33) 



The pilot signals are made up of segments, each having a duration T s , 

referred to as the pilot symbol duration, over which the pilot signals are 
10 orthogonal, and wherein the following holds 



15 



20 



25 



J x* (t) ■ Xj (t)dt = 0 V i, j 6 {IK N Tx } i*j 

(n-\)T s 



(34) 



The pilot vector is defined as 



2(f) = 



x,(t) 



(35) 



and the receiver noise signal z(t) represents white Gaussian noise. The signal 



received by the single Rx antenna is described as 

y(t) = 8 T <f).K(t) + z(t)., 



(36) 



Conventionally, correlating the received signal with the four (4) pilot 
sequences derives a set of four (4) channel estimates. Wherein the pilot signals 
are orthogonal over a pilot symbol period, this estimation is then repeated at 
the pilot symbol rate. Such a correlation procedure is generally referred to as 
"integrate and dump" and may be expressed as 



(37) 



S (n-l)T s 
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wherein h™ v is a vector made up of conventional, i.e., integrate and dump, 

channel estimates derived from the n-th pilot symbol. If (34) is transformed 
into a discrete time representation, by putting N T =T S IT samples of the pilot 

signals x n (t) into the columns of the matrix A , N T samples of the noise signal 

5 z(t) into the vector K , and N T samples of the received signal y[t) into the 
vector ¥ , (34) yields 

?= A-fe> + #. (38) 



10 



Then the discrete time representation of (37) is 



Jfw = J_A H -P. (39) 
N T 



cov 



If the channel variations within one pilot symbol are neglected, (39) becomes 

15 fe=£(nr s ) + -L.A H .£. (40) 

Considering the linear transformation of h(t) into %{t), the received 
signal is expressed as 

20 ¥=A%{a)h{t)+H. (41) 

From this, an estimate of the scalar h(t) is derived as 



= SVpX m Jfw (42) 

c ° nv 

25 

Again, when the channel variations within one pilot symbol are ignored, (13) 
becomes 

h"-h(«T s ) + 8 * T J*\ ■ A" • ft '. (43) 
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From this scalar estimate, using the linear transformation, a new estimate of the 
channel impulse vector is generated as 



fe-ifr).** = %^; (a) -fe,- (44) 



Ignoring channel variations within one pilot symbol, (44) becomes 



&=f(,r l)+ teg.A".f. (45) 



10 If £{cc) is not known a-priori, it may be estimated using the covariance 

matrix given by 



R h = (K(nT s ) ■ K* (nT s )) = 8(a) - («)•/>, (46) 



15 with F A being the average power of the scalar channel impulse response hit) 
The covariance matrix R h can be approximated as 



^=77-Zfcv-fcv H , (47) 



^ sym «=1 



20 which averages the vector with the conventional channel impulse response 
estimates over a number, N sym , of pilot symbols. 

For the case without noise and having an angular spread equal to zero, 

R h is rank one (1) and the vector $(a) spans R h . Thus (47) reduces to 



25 R 



h 



noiseless 



= J(a)^(a) 



Note that the normalized vector &{<z)/\\ii(czj\ spans R a . 

For a noisy case with sufficient low noise power and sufficient low 
angular spread, is still dominated by one eigenvalue. Therefore, the process 

30 performs an eigenvalue decomposition of R h . When one eigenvalue is much 
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larger than all other eigenvalues, it is an indication that the angular spread 
around was rather small. Therefore, as f max is the eigenvector 

corresponding to the largest eigenvalue of R h , the approximation becomes 

fa) - g <r (a) $ $*r (49) 

a II P/ Ml 2 max ^max * V*- 7 J 

Note that the vector in this example is equal to the channel subspace E c . 
In general, the estimate R h is used to determine whether the rank of the channel 

estimation covariance matrix can be reduced. If R h is full rank, the channel 

10 estimation problem is not reduced to a smaller dimension. 

According to the exemplary embodiment, orthogonal pilot signals of 
binary chips have a chip rate of 1.2288 Mcps, and a pilot symbol duration of 64 
chips. With this channel model, a received signal, including white Gaussian 
noise, is generated for 4000 consecutive pilot symbols having a pilot Signal-to- 

15 Noise Ratio (SNR). From the received signal, 4000 conventional vector 

estimates, fi™ v , are generated. The thus generated covariance matrix R h is 

averaged over these 4000 consecutive conventional channel estimates. In the 
exemplary embodiment, the process takes approximately 208.3 ms, After 

extracting the eigenvector corresponding to the maximal eigenvalue of R h , the 
20 matrix R a is calculated. Subsequently, 4000 new vector estimates ft™ are 



produced according to 



fc = Ra ' ki • (50) 



25 Using the exemplary embodiment, iterations are repeated N exp = 50 

times. Over the 50 iterations the transmitter angles are varied such that a is 
uniformly distributed within (+/- 60) degrees, while the angular spread 

remains constant, having a standard deviation of square root of two (V2) 
degrees. Additionally, the channel parameters for a given pilot SNR are varied. 



WO 02/060138 



PCT/US01/50068 



26 

The varied parameters represent radio path direction(s), path-specific phase, 
and path-specific Doppler shift, for a certain pilot SNR. An equal number of 
iterations is performed for different pilot SNR values. A comparison of the 
quality of the set of conventional estimates to the set of new vector estimates, 
5 with respect to the reduction factor of the mean squared estimation error that is 
averaged over time and iterations, is made using the estimation gain given as 




10 For the exemplary embodiment, FIG. 6 illustrates the estimation gain in 

dB as a function of the pilot SNR. Wherein the pilot SNR is defined as the ratio 
of the average energy per pilot chip E c of one pilot signal received at the single 

antenna receiver to the received noise power density I 0 in dB. 

The upper limit for the estimation gain is determined by the number of 
15 transmit antennas, which is illustrated in FIG. 6 as 6dB. As illustrated in FIG. 6, 
the estimation gain approaches the upper limit even though the assumed 
angular spread is not zero and the received signal is severely corrupted by 
noise. The reduction of the estimation gain with increasing pilot SNR is due to 
the non-zero angular spread. 
20 Although the channel impulse responses are not completely correlated, 

the derivation of the impulse response n new assumes this property. For larger 

angular spreads, a smaller estimation gain is expected. For small angular 
spreads, the estimation gain appears considerable. Note that in general, for 
residential and suburban environments a standard deviation of one (1) to two 
25 (2) degrees is frequently observed. Note also that it is possible to evaluate the 
performance improvement of the reduced rank channel estimation method 
using a Monte-Carlo-simulation to derive the reduction of estimation errors of 
the channel impulse responses as compared to conventional channel estimation 
using independent correlators. 
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Reduced rank channel estimation for systems using multiple transmitter 
antennas allows improvement of the channel estimation quality under certain 
propagation conditions with limited diversity due to correlated fading. As the 
mechanisms affecting correlation, such as the directionality of the radio wave 
5 propagation, change relatively slowly over time, the correlation characteristics 
may be estimated by averaging over extended time intervals. This is in contrast 
to the time intervals associated with inverse fading rate of the channel and thus 
allows improved accuracy in estimating the correlation characteristics. 

Reduced rank channel estimation for multiple transmitter antennas is 
10 also applicable to frequency-selective channels by computing either separate 
estimates of the correlation characteristics or by computing estimates of the 
correlation characteristics across all propagation delays. Separate estimates 

refers to computation of R h , for each propagation delay. Reduced rank 

channel estimation is then performed taking into account each delay occurring 
15 in the frequency-selective channel impulse response. In an alternate 
embodiment, wherein additional information, such as the antenna 
configuration at the transmitter, is known a-priori, the step of estimating the 
linear transformation of the reduced number of uncorrelated channels into the 
larger number of correlated channels may be more accurate. Additionally, the 
20 reduced rank estimation process may be extended to cases with more than one 
receiver antenna. In this case, the estimation is performed for the MIMO 
channels, as illustrated in FIG. 1. While the present example involves a system 
employing coherent demodulation, reduced rank channel estimation as 
described herein is also applicable to communication systems employing non- 
25 coherent demodulation. 

A receiver 100 according to one embodiment of the present invention is 
illustrated in FIG. 7. The receiver 100 has a single antenna 102 that receives 
signals from a transmitter having multiple antennas. The received signals are 
first processed by the preprocessor 104. The signals are then provided to a 
30 correlator 106, which is used as a sliding correlator for searching and as a 
correlator for the significant delays for channel estimation. In an alternate 



WO 02/060138 



PCT/USO 1/50068 



28 

embodiment the delays are determined in software without use of a correlator. 
The outputs of the correlator 106 are used to provide an estimate of the 
covariance matrix. In one embodiment, the correlator 106 is made up of fingers 
to form a rake, having one finger for each combination of transmitter antenna, 
5 receiver antenna and significant delay. The estimates are provided to the 
central processor 112 via bus 116. The processor 112 stores the channel 
parameter estimates in memory 114 so that the estimates may be used to derive 
the channel covariance matrix averaged over time. 

From memory 114, the estimated covariance matrix is provided to the 

10 rank analysis and subspace estimation unit 108 for eigenvalue decomposition. 
If one or more eigenvalues dominate the others, the channel subspace is 
estimated by computing the eigenvectors that correspond to the dominant 
eigenvalues. The eigenvectors spanning the channel subspace are written to 
memory for further use in the channel subspace projection unit 109 where 

15 reduced rank channel parameter estimates are produced by computing the 
projection of the (N Tx • • N E ) original channel estimates per estimation time 

interval onto the channel subspace, yielding N Ch reduced rank channel 

parameter estimates per estimation time interval. The results of the channel 
subspace projection unit 109 are written to memory for use in the demodulator 
20 110. Optionally the channel subspace projection unit 109 could generate 
equivalent full dimension channel parameter estimates, by re-transforming the 
N Ch reduced rank channel parameter estimates into (N Tx • -N E ) equivalent 

full dimension channel parameter estimates per estimation time interval. For 
example in a conventional RAKE-receiver design for the full rank channel 
25 model, the number of rake fingers for a full rank demodulator would be 
(N Tx • • N E ) . A full rank demodulator would then use the (N Tx • • N E ) 
original channel parameter estimates for the finger coefficients. A reduced 
complexity demodulator could eventually use only N Ch RAKE fingers using the 

N Ch reduced rank channel estimates as coefficients. However, since the 
30 receiver would generally be designed in anticipation of a worst case situation, 
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i.e., wherein (N Tx • -N E ) fingers are implemented, it would be sufficient to 
compute ( N Tx • * N E ) correlated channel parameter estimates with 
improved estimation quality over of the N ch reduced rank channel parameter 
estimates. 

5 The rank analysis and subspace estimation unit 108 and the subspace 

projection unit 109 may be implemented in a Digital Signal Processor (DSP), 
dedicated hardware, software, firmware, or a combination thereof. Modules 
within receiver 100 may be incorporated together, and are illustrated as 
separate blocks for clarity based on function. 

10 An exemplary configuration of one embodiment is illustrated in FIG. 8 

for a system having four (4) transmitter antennas and two (2) receiver antennas. 
Three (3) transmission paths are illustrated and labeled 1, 2 and 3. The points of 
reflection for paths 1 and 2 are both on a same ellipse, wherein the ellipse is 
formed such that Tx and Rx are the focal points. Note that the ellipse is 

15 superimposed on the illustration of the physical layout of the system. Path 3 
falls outside of the illustrated ellipse. Paths 1 and 2 have the same significant 
delay, Xi, with respect to the receiver, while path 3 has a significant delay, x 2 
different from x,. The path delay is a function of the configuration of the 
antennas as well as the environment of the system. As illustrated, the four (4) 

20 transmitter antennas and the two (2) receiver antennas result in eight (8) 
channels. Each of the path delays, i\ and x 2 , produce an echo, wherein ( N E =2). 
The dimension of the covariance matrix is given as ( N Tx • • N E ) or sixteen 
(16) corresponding to the ( N Tx • * N E ) channel impulse response samples. 

Therefore, the full rank channel parameter vector is a 16-dimension vector. 
25 Using the rank reduction methods described herein, the rank of the channel 
estimation may be reduced to three (3) dimensions, corresponding to paths 1, 2, 
and 3, wherein (A^ =2). Note that where the mapping of the N Ch 

transmission paths to the ( N Tx • • N E ) channel impulse response samples is 

not known, the subspace may be extracted from configuration information. If 
30 the location and characteristics, such as direction and directionality, of the 
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antennas are known, the information may be used to generate an array 
response or steering vector. Using the steering vector and path direction 
information, which is also extractable using subspace algorithms, the angle of 
transmission, a, is estimated. If the antenna configuration has a fixed 
5 deployment the angle of transmission is calculable. A vector is formed 
including an angle of transmission for each transmitter antenna. Similarly, an 
arrival angle vector is formed considering the receiver antennas. A linear 
transformation for the mapping of the N Ch transmission paths to the 

( N Tx - N to - N E ) channel impulse response samples is constructed using this 

10 information from both the transmitter and receiver configurations. This 
provides the matrix B as given in (19) hereinabove describing the linear 
transformation. The covariance matrix is derived therefrom as in (20) 
hereinabove. The process then proceeds as for the case where corresponding 
information is obtained from a priori knowledge. 

15 While one embodiment has been described herein with respect to the 

time domain, an alternate embodiment performs a rank reduction of the 
covariance matrix or a sample matrix in the frequency domain. If the 
parameters and equations are developed in the frequency domain, the process 
to estimate the channel then incorporates the frequency domain values. 

20 The previous description of the preferred embodiments is provided to 

enable any person skilled in the art to make or use the present invention. The 
various modifications to these embodiments will be readily apparent to those 
skilled in the art, and the generic principles defined herein may be applied to 
other embodiments without the use of the inventive faculty. Thus, the present 

25 invention is not intended to be limited to the embodiments shown herein but is 
to be accorded the widest scope consistent with the principles and novel 
features disclosed herein. 
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CLAIMS 

I (WE) CLAIM: 

1. A method for modeling a link in a wireless communication system, the 
2 system having a transmitter having N antennas and a receiver having M 

antennas, each path from one of the N transmitter antennas to the M receiver 
4 antennas comprising a channel, the method comprising: 

determining a matrix describing parametric relations of the link; 
6 ranking the matrix; 

determining if the rank is less than NxM; 
8 if the rank is less then NxM performing an extraction of a subspace of 

the matrix; 

10 deriving channel impulse responses for each channel based on the 

extracted subspace of the matrix; and 
12 demodulating a received signal using the channel impulse responses. 

2. The method of claim 1, wherein the matrix is a covariance matrix 

2 describing the link, wherein the covariance matrix represents a plurality of 
impulse responses between the transmitter and the receiver. 

3. The method of claim 1, wherein the matrix is a sample matrix describing 
2 the link. 

4. The method of claim 1, wherein the step of determining the matrix 
2 further comprises: 

estimating a plurality of parameters describing at least one channel. 

5. The method of claim 4, wherein the parameters include a distance 
2 between transmitter antennas. 
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2 6. The method of claim 4, wherein the parameters include a transmittal 
angle with respect to a configuration of the transmitter antennas. 

7. The method of claim 4, wherein the determining the matrix comprises 
2 estimating the matrix. 

8. The method of claim 1, wherein the matrix describes parametric 
2 relations of the link in the frequency domain. 

9. The method of claim 1, wherein the ranking the matrix, further 
2 comprises: 

determining an eigenvalue for the matrix. 

10. The method of claim 1, wherein if the rank is equal to (NxM) a set of 
2 correlated impulse responses is applied for demodulating. 

11. A wireless apparatus operative to perform the method of claim 1. 

12. A wireless communication apparatus, comprising: 

2 a correlator operative to estimate a covariance matrix representing a link 

with a transmitter based on signals received from the transmitter; 

4 a rank analysis unit coupled to the correlator and operative to estimate a 

rank of the covariance matrix; and 

6 a channel estimation unit coupled to the rank analysis unit and operative 

to generate a reduced rank channel estimate. 

13. The apparatus of claim 12, wherein the covariance matrix represents a 
2 plurality of impulse responses between the apparatus and the transmitter. 

14. The apparatus of claim 12, wherein correlator is operative to determine a 
2 correlation of at least two channels. 
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15. The apparatus of claim 14, wherein the rank analysis unit is operative to 
2 determine an eigenvalue corresponding to the co variance matrix. 

16. The apparatus of claim 15, wherein the rank analysis unit is operative to 
2 compare the estimated rank of the covariance matrix to a predetermined full 

value. 

17. A method for estimating a link in a wireless communication system, the 
2 method comprising: 

estimating a covariance matrix for the link; 
4 determining if the rank of the covariance matrix is reducible; 

reducing the rank of the covariance matrix; and 
6 estimating a set of impulse responses for the link using the reduced rank 

covariance matrix. 

18. The method of claim 17, further comprising: 
2 determining a correlation of the channel; 

ranking the covariance matrix; and 
4 performing an extraction of a reduced rank matrix out of the covariance 

matrix. 

19. A wireless communication apparatus operative within a wireless 

2 communication system having a transmitter having N antennas and a receiver 
having M antennas, each path from one of the N transmitter antennas to the M 
4 receiver antennas comprising a channel, the apparatus comprising: 

a first set of computer readable instructions operative to determine a 
6 covariance matrix describing the link; 

a second set of computer readable instructions operative to rank the 
8 covariance matrix; 

a third set of computer readable instructions operative to determine if 
10 the rank is less than NxM; 
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a fourth set of computer readable instructions operative to perform an 
12 extraction of a reduced rank matrix out of the covariance matrix 

if the rank is less then NxM; 
14 a fifth set of computer readable instructions operative to derive channel 

impulse responses for each channel based on the reduced rank 
16 covariance matrix; 

a sixth set of computer readable instructions operative to demodulate a 
18 received signal using the channel impulse responses. 



20. The apparatus of claim 19, further comprising: 

an equalizer operative in response to the sixth set of computer readable 
instructions, wherein a configuration of the equalizer is 
determined by the rank of the covariance matrix. 



21 . The apparatus of claim 19, further comprising: 
2 a seventh set of computer readable instructions operative to derive a correlated 
channel impulse response. 



22. A wireless communication apparatus, comprising: 
2 a channel estimation means operative to estimate a covariance matrix 

representing a link with a transmitter based on signals received 
4 from the transmitter; 

a rank analysis unit coupled to the correlator and operative to estimate 
6 the rank of the covariance matrix; and 

a channel estimation means coupled to the rank analysis unit and 
8 operative to generate a reduced rank channel estimate. 



23. A wireless communication apparatus, comprising: 

a correlator operative to estimate a covariance matrix representing a link 
with a transmitter based on signals received from the transmitter; 
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4 a rank analysis unit coupled to the correlator and operative to estimate 

the rank of the co variance matrix; and 
6 a channel estimation means coupled to the rank analysis unit and 

operative to generate a reduced rank channel estimate. 

24. A method for estimating a link in a wireless communication system, the 
2 method comprising: 

estimating a covariance matrix for the link; 
4 determining if the rank of the covariance matrix is reducible; 

reducing the rank of the covariance matrix; and 
6 estimating a set of impulse responses for the link using the reduced rank 

covariance matrix. 



25. The method of claim 24, further comprising: 
2 determining a correlation of the channel; 

ranking the covariance matrix; and 
4 performing an extraction of a reduced rank matrix out of the covariance 

matrix; 



26. A wireless apparatus, comprising: 
2 channel estimation means operative to determine significant delays and 

determine a set of estimates of full dimension channel parameters 
4 associated with the significant delays, wherein each one of the set 

of estimates corresponds to an instance in time; 
6 eigenvalue computation means operative to determine eigenvalues of 

the set of estimates of the full dimension channel parameters and 
8 find any dominant eigenvalues; and 

channel estimation means operative to determine a set of reduced rank 
10 channel parameter estimates in response to the dominant 

eigenvalues. 
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27. The wireless apparatus of claim 26, further comprising: 

eigenvector computation means operative to determine at least one 
eigenvector associated with one of the dominant eigenvalues of 
the set of estimates; 
wherein the channel estimation means uses the at least one eigenvector 
to project the set of estimates of the full dimension channel 
parameters onto the subspace spanned by the at least one 
eigenvector. 
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